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Abstract. There has been a recent stimulus in the study of alternative theories of gravity 
lately, mostly triggered from combined motivation coming from cosmology /astrophysics and 
high energy physics. Among the proposed theories, one that has attracted much attention 
, is f{R) gravity. It is certainly debatable whether such a simplistic modification of General 

' Relativity can constitute a viable alternative theory of gravitation. However, it is quite 

straightforward to see the merits of such a theory when viewed as a toy theory whose role 
is to help us understand the implications and difficulties of beyond-Einstein gravity. Under this 
perspective, I review some of the main lessons we seem to have learned from the study of f{R) 
gravity in the recent past. 



^ I 1. Introduction 

' Modifications of General Relativity (GR) by including higher order curvature invariants in 

l/^ . the gravitational action have a long history. Indeed, it was just 1919 when Weyl, and 1922 

I when Eddington (the very man that three years earlier had provided the first experimental 

verification of GR by measuring light bending during a solar eclipse) started considering such 
terms in the action [H [2] . These early attempts were mostly driven by scientific curiosity and a 
need to question and, therefore, understand the then newly proposed theory. Thus, they were 
quickly abandoned. However, by the 1960's and 1970's new motivation coming from theoretical 
physics revived the field of higher-order gravity theories: Utiyama and De Witt showed that 
^ ■ renormalization at one-loop demands that the Einstein-Hilbert action be supplemented by higher 

order curvature terms p]. Later on, Stelle showed that higher order actions are renormalizable 
(but not unitary) |1] (see [5] for a historical review and a list of references to early work) . More 
recent results show that when quantum corrections or string theory are taken into account, 
the effective low energy gravitational action admits higher order curvature invariants [SI [71 18] . 
However, the relevance of such terms in the action was considered to be restricted to very strong 
gravity regimes and they were expected to be strongly suppressed by small couplings, as one 
would expect when simple effective field theory considerations are taken into account. 

Recently, a new stimulus appeared in higher-order theories of gravity. This time the 
motivation came from cosmological and astrophysical observations. The latest datasets coming 
from different sources, such as the Cosmic Microwave Background Radiation (CMBR) and 
supernovae surveys, seem to indicate that the energy budget of the universe is 4% ordinary 
baryonic matter, 20% dark matter and 76% dark energy [9[ \10 \ [TT | I12j. The term dark matter 
refers to an unkown form of matter, which has the clustering properties of ordinary matter but 
has not yet been detected in the laboratory. The term dark energy is reserved for an unknown 
form of energy which not only has not been detected directly, but also does not cluster as ordinary 



matter does. The simplest model that fits the data adequately, is the so called concordance model 
or ACDM model (supplemented with some inflationary scenario) , where a cosmological constant 
plays the role of dark energy and the nature of dark matter remains undetermined apart from 
the fact that it does not interact but gravitationally. However, even in the most optimistic 
scenario in which the nature of dark matter will soon be determined in some accelerator, the 
known cosmological constant problems [13^ [H] would still plague the concordance model. One 
of the alternatives is to consider the unexpected observation as shortcoming of our macroscopic 
description of gravity itself: gravity should be by far the dominant interaction at cosmological 
scales so, could it be that a modification of GR could account for the unexplained observations? 

It is tempting to combine the earlier motivation for corrections to GR coming from high 
energy physics with the more recent cosmological/astrophysical stimulus. This is indeed often 
done in the literature by proposing that some theory of gravity, which comes about as a low 
energy limit of a more fundamental theory and includes both ultraviolet and infrared corrections, 
could maybe account for all or some of the unexplained observations. However, there are two 
loopholes in such a suggestion: a) It is not clear exactly which one the fundamental theory 
should be or what is its low energy limit, and b) the length scales at which we expect high- 
energy inspired corrections to GR to be important are orders of magnitude smaller than the 
length scales at which our puzzling observations have been made. 

These two issues can be seen as caveats against modified gravity; or, in a more optimistic 
approach, they can be seen as opportunities: Problems, such as the "unnatural" value of the 
cosmological constant have not found a solution within the realms of conventional (or even 
slightly unconventional) thinking. Therefore, they might as well be indications that something 
is wrong with the very core of this way of thinking. The large ambiguities on how gravity 
works at small scales and which is the fundamental theory describing it seem to leave enough 
room for such speculations (indeed there are results already in the literature claiming that 
terms responsible for late time gravitational phenomenology might be predicted by some more 
fundamental theory, such as string theory [15] ) . 

What seems to be the best way to proceed is probably the usual way when we are diving into 
the unknown in physics: consider a toy theory, probably of a very simple form, and use it as a 
tool to explore the limitations of the theory being questioned, in our case GR. This would allow 
us to test whether the very idea of modified gravity accounting for unexplained cosmological 
observations can be correct in principle, or what is the relation between the modifications 
required for that and those predicted by high energy physics. 

/(i?) theories of gravity seem to be very good candidates for such toy theories. They come 
about by straightforwardly generalizing the Lagrangian in the Einstein-Hilbert lagrangian to be 
a general function / of the Ricci scalar R, so that the action becomes 



where k = SirG, G is the gravitational constant, g is the determinant of the metric and R is 
the Ricci scalar [c = h = 1). f{R) theories of gravity make good toy theories for two reasons: 
a) They are sufficiently general to encapsulate some of the basic characteristics of higher-order 
gravity, but at the same time they are simple enough to be easy to handle; b) there are unique 
among higher-order gravity theories, in the sense that they seem to be the only ones which can 
avoid the long known and fatal Ostrogradski instability [16j. 

Now, if we are to consider f{R) gravity as a toy theory whose main scope was to help us 
learn something about beyond GR gravity, the obvious question to ask is what did we really 
learn so far? I attempt to answer this question here by going through some "lessons" coming 
from f{R) gravity (hence, the title). This talk is based on [T7], where a longer discussion on the 
motivation of f{R) gravity, as well as all on the other issues that are discussed here, is given 
(see also [18]). 





2. What have we learned from f{R) gravity? 

2.1. Lesson 1: Choice of fundamental variables is crucial 

There are actually three versions of f{R) gravity: Metric f{R) gravity, Palatini f{R) gravity, 
and metric-affine f{R) gravity. In metric f{R) gravity, the action 

Sraet = ^ j d^^^^-Q fi^) + SM{9fiu,1p), (2) 

where a matter action Sm has been included and ip collectively denotes the matter fields, is 
varied with respect to the metric as usual. This yields the field equations 

f'(.R)R^.. - \f{R)gp.u - [V^V, - g^M f'{R) = kT^,, (3) 

where 

a prime denotes differentiation with respect to the argument, is the covariant derivative 
associated with the Levi-Civita connection of the metric, and □ = V^V^. Clearly, these equation 
are fourth order differential equations in the metric. 

Palatini f{R) gravity comes about from the same action if we decide to treat the connection 
as an independent quantity. In this case, for the sake of clarity, we denote as TZ^i, the Ricci 
tensor constructed with this independent connection. TZ = g^'^TZ^j^u is the corresponding Ricci 
scalar and the action takes the form: 

^P<^i ^ ^ j ^'^^"'^ f^'^'> ^ ^^''^3iiu,tp)- (5) 

Independent variation with respect to the metric and the connection gives, after some 
manipulations, the following field equations 

/'(7^)7^(,.,) - ^/(7^)5;.. = kGv, (6) 

Va (^/'(7^)5^-) = 0, (7) 

where T^i, is defined in the usual way as in eq. (jH) and denotes the covariant derivative 
defined with the independent connection. Notice that when f{TZ) = TZ, eq. d?]) becomes the 
definition of the Levi-Civita connection and eq. © becomes Einstein's equation. So Palatini 
and metric variations both lead to GR for an action linear in the Ricci scalar. However, they 
lead to different theories for a more general action. In this sense, it is not straightforward which 
variation on should choose when generalizing the lagrangian to be a general function of the Ricci 
scalar. Both metric and Palatini f{R) gravity was first rigorously studied in |19j. 

In Palatini f{R) gravity the connection, although independent, does not enter the matter 
action (this means that Palatini f{R) gravity is a metric theory, according to the definition of 
[20j ) . This implies that the covariant derivative of the matter fields which are usually present in 
the matter action are defined with the Levi-Civita connection (alternatively, the theory would 
include only very specific matter which does not couple to the connection, such as scalar fields, 
but this is not really an option for a gravity theory) [211 [22l [23] . Additionally, the independent 
connection is really an auxiliary field: eq. ([7]) can be solved algebraically for the connection to 
give 



d,,{nn)g,,)+d,{f'{n)g 



da{f'{n)g^,) . (8) 



At the same time, a contraction of eq. ([6]) with the metric gives 

nn)n-2f{n) = KT, (9) 

where T = g^'^Tnj,. for a given / this is an algebraic equation which can be solved to express 
TZ in terms of TjH Thus, using eqs. dHj) and ([9]) one can algebraically express the independent 
connection in terms of the metric and the matter fields and eliminate it from the field equations. 

The way to give more substance to the connection is to restore its geometrical role and use 
it to define the covariant derivatives in the matter action, as proposed in [22j. In this case the 
matter action would be of the form 5*/ (g'^jy, F"^^^, ip). Such theories are called metric-afhne f{R) 
theories. We refrain from giving their field equations here as we will not use them in what comes 
next. We refer the reader to [22] for more details. 

What should be clear after this discussion is that it is crucial to pick, not only the action, 
but also the fundamental fields which will be varied and assign to them a precise meaning. 
Remarkably, the distinction between Palatini and metric variation, which (at least at a classical 
level) was of little importance for the Einstein-Hilbert action, is critical for more general actions. 



2.2. Lesson 2: Equivalence with Brans-Dicke theory 

Metric f{R) gravity appears as an extension of GR that includes no extra fields mediating 
gravity. Palatini f{R) gravity appears to include a independent connection, so one could be 
trick to think that the underlying geometry is not pseudo-Riemannian. We already argued 
that the independent connection is not really related to the geometry (at least as the latter 
is felt by the matter field which act as our probe). So appearances can deceive; and they 
do: both metric and Palatini f{R) can be written as specific versions of Brans-Dicke theory 
[251 ESI EZl [Ml Eg [301 [SI ED El 123] . Let us see how this is possible. 

Starting from the action ([2|) for metric f{R) gravity, one can introduce a new field x ^'^^ 
write the dynamically equivalent action 

Smet = Y^j d'xV^ [fix) + f'{x){R - x)] SM{g^u,^)- (10) 

Variation with respect to x leads to the equation 

r(x)(i?-x) = o. (11) 

Therefore, x = ii if f"{x) 0) which reproduces the action ([2|). Redefining the field x by 
'P = fix) and setting 

Vicl>) = xm-fixm, (12) 

the action takes the form 

Smet = ^ / d^^V^ [<PR - Vm + SM{g^^., V')- (13) 

This is the Jordan frame representation of the action of a Brans-Dicke theory with Brans-Dicke 
parameter ujq = (also called "O'Hanlon action" as it was originally proposed by O'Hanlon in 
[33]). 

We could follow exactly the same steps starting from action dSj) for Palatini f{R) gravity. We 
would then arrive to the action 

Spal = ^ I d^X^ [<pn - Vm + SM{9^u,i^). (14) 



^ See |24l 117] for a discussion on the exceptional case where theis equation has no solutions or when f(TZ) oc TZ^ . 



This action is identical to (I13p apart from the fact that R has been now replaced by TZ. This 
is critical, as TZ is not the Ricci scalar of the metric g^u and, therefore, action (I14p is not that 
of a Brans-Dicke theory with loq = 0. However, taking into account also the definition of (p, we 
known that we can express the independent connection algebraically with respect to the metric 
and (p using eq. El Using this expression action ([Tl]) takes the form 

Spai = Y^j (fx^g (^^R + ^d^<t^d^^ - V{<P)j + SMig^u, (15) 

This is the action of a Brans-Dicke theory with Brans-Dicke parameter = —3/2. 

I will avoid to present the complete set of field equations for Brans-Dicke theory for the sake 
of brevity and give only the equation governing the dynamics of the scalar: 

{2u)Q + 2,)U(t) + 2V - (pV = kT. (16) 

This equation is enough to reveal that, it is not the ujq = case the one for which (j) carries no 
dynamics as one might have expected judging from the absence of a kinetic term in the action, 
but actually the cjq = —3/2, which corresponds to Palatini f{R) gravity. 

The outcome of this discussion is twofold: Firstly, both metric and Palatini f{R) gravity 
are equivalent to some version of Brans-Dicke theory (with all the implication this carries: for 
instance, this highlights the fact that Palatini f{R) gravity is really a metric theory). Secondly, 
in metric f{R) gravity there is only one extra scalar degree of freedom, where as in Palatini, 
the Brans-Dicke scalar is non-dynamical (so it is an auxiliary field, much like the independent 
connection). The moral is that not all representations of a theory are equally suitable for 
straightforwardly exhibiting some of the features of the theory (for a longer discussion see also 
iMj). 



2.3. Lesson 3: Devil hides in the details 

During the last few years it has been understood that it is not that difficult to construct models 
within the framework of f{R) gravity that lead to a desired cosmological background evolution 
(see [MllMlEIllMllMlllOlSIlllSlllSllHlligfor some examples). Without getting into the details, 
the procedure is quite simple: Under the usual cosmological assumptions regarding symmetry 
and matter description, f{R) gravity leads to modified Friedmann equations governing the 
dynamics of the universe. Choosing the function / properly can lead to the desirable modification 
of the dynamics {e.g. late time expansion). In fact, it has even been shown that, starting from 
a given scale factor a{t) as a function of cosmic time, one could "reconstruct" the function / (or 
at least a family of functions /) that leads to this evolution }46l Wl\ 148] . 

However, the difficulties arise when one goes beyond background evolution and requires that 
some f{R) can account also for the fine details in cosmological observations. For instance, the 
requirement that de Sitter spacetime (in models which do not admit Minkowski space as a 
maximally symmetric solution) be a stable solution gives important constraints [49ll50l[5T] . Big 
Bang Nucleosynthesis is a delicate process that can lead to severe constraints as well {e.g. [52j). 
The growth of cosmological perturbation and structure formation are other examples of processes 
that can rule out models and lead to strong viability constraints [ 531 [Slj [55] . See p!7] for a more 
detailed discussion on cosmological constraints and a more complete list of references. 

The lesson learned here is that the devil always hides in the details and proposing 
cosmologically viable f{R) gravity models can be significantly trickier than it seems. The subtle 
details of the cosmological evolution are quite sensitive to the underlying gravity theory and 
as we are entering the era of precision cosmology it should become more and more difficult to 
deviate from GR without affecting the delicate picture of the standard model of cosmology^] 



^ Besides cosmology, another example where perturbations around the spacetime are much more revealing than 



2.4- Lesson 4-' Extra degrees of freedom are troublesome 

As mentioned in the Introduction, metric f{R) gravity can avoid the Ostrogradski instability, 
unUke higher order gravity theories in general |16j . At the same time f{R) theories have no 
ghosts (massive states of negative norm that cause apparent lack of unitarity), which is also not 
true for generic higher order theories of gravity [U El [58l El HI [Ml EQ]. Therefore, one could say 
that metric f{R) gravity has been carefully chosen so that the extra degrees of freedom with 
respect to those of GR will not evidently lead to serious viability issues. 

Nevertheless, these criteria have proven inadequate and problems related with the presence 
of extra degrees of freedom are still plaguing metric f{R) gravity. The most prominent is the 
instability initially found by Dolgov and Kawasaki [HI] for the model proposed in |36j but can 
appear for more general functions / [62j|f| To see how this comes about let us parametrize the 
deviations from Einstein gravity as 

f{R) = R + ev{R), (17) 

where e is a small parameter with the dimensions of a mass squared and (p is arranged to be 
dimensionless. The trace of eq. dSj) yields 

3af'{R) + f'{R)R-2f{R) = KT, (18) 

and substituting eq. (fT?!) gives 



□fl+^V<.;;V„fl+<if^ii=^ + |^. (19) 

ip" Se if ' ie if ' 6ip 

where it is assumed that ip" ^ 0. Consider now a small region of spacetime in the weak- field 
regime and approximate locally the metric and the curvature by 

gfiu = ViMu + h^u, R = -kT + Ri, (20) 

where rj^n is the Minkowski metric and \Ri/kT\ <^ 1. Equation (f20l) yields, to first order in Ri, 

Ri - V^Ri ^ TRi + — ^ VT • VRi 

ip" ip" 

1 A An rj^ V72^ {KTip' + 2ip) 



-(-- ip'] Ri = kT-kV^T 
o" \e ) 



Zip" \e ) 2,ip" 

(21) 

where V and are the gradient and Laplacian in Euclidean three-dimensional space, 
respectively, and an overdot denotes differentiation with respect to time. The function ip and 
its derivatives are now evaluated at = —kT. The coefficient of Ri in the fifth term on the 
left hand side is the square of an effective mass and is dominated by the term (3e(/3")~^ due to 
the extremely small value of e needed for these theories to reproduce the correct cosmological 
dynamics. Then, the scalar mode Ri of the f{R) theory is stable if ip" = f" > 0, and unstable 
if this effective mass is negative, i.e., if ip" = f" < 0. The time scale for this instability to 
manifest is estimated to be of the order of the inverse effective mass ~ 10~^^ s in the example 

the background itself, is black holes: even though the black hole solutions of GR are also solutions of f{R) 
gravity model and one cannot distinguish between theories by mapping the spacetime around a black hole [56], 
the perturbations around a black hole solution behave differently in different theories and, therefore, they allow 
for distinctions [57] . 

This instability does not appear in Palatini f(R) gravity [63] , 



eip{R) = —iJL^/R [61j- The small value of ip" gives a large effective mass and is responsible for 
the small time scale over which the instability develops. 

This fatal instability for models with /" < manifests itself in the linearized version of 
equation ()18p . There are also recent claims that R can be driven to infinity due to strong non- 
linear effects related to the same equation [641 [65l [66] . Both of these problems are strictly related 
to the extra scalar degree of freedom in metric f{R) gravity and hint towards the following 
conclusion: no matter how carefully one constructs a theory with extra degrees of freedom, 
complications related to them can arise as extra degrees of freedom are very difficult to control. 

2. 5. Lesson 5: Parametrized post- Newtonian expansion is difficult hut fruitful 
Even though one might think that the parametrized post-Newtonian (PPN) expansion is 
more or less a standardized procedure, it actually turned out to hold a lot of surprises 
in f{R) gravity. It took several papers and lengthy debates for this issue to be settled 
{e.g. [SZllMlESlEIllTOlEIlEaEa], see also ^ for a more detailed hst of references). The 
main problems which were encountered and led to debates were: a) the equivalence with Brans- 
Dicke theory and whether this can be used for the PPN expansion, b) confusion related to the 
absence of flat spacetime as a solution in many models, c) misunderstandings regarding the 
difference between the existence of the Schwarzschild(-de Sitter) metric as a vacuum solution 
and a correct (post-)Newtonian limit, d) difficulty with non-analytic functions f{R) and their 
expansion. Most of these problems were not specific to f{R) gravity and, therefore, solving them 
has led to a better understanding of the post-Newtonian expansion for alternative theories of 
gravity in general. 

In spite of its difficulties, the post-Newtonian limit of f{R) gravity has proved to be very 
fruitful after all: it led to strong viability constraints. Briefly summarized, the outcome of the 
relevant studies is the following: For metric f{R) gravity, one can indeed use the equivalence 
with Brans-Dicke theory, but with some care {e.g. the scalar is algebraically related to the 
curvature so it does not tend to settle to the minimum of the potential). This implies that, since 
the theory corresponds to a Brans-Dicke parameter ujq = 0, it is ruled out (experiment require 
Wo > 40000) unless the effective mass of the scalar is large. The effective mass is essentially the 
second derivative of the potential whose functional form depends on the form of /. It has also 
been shown that certain models of metric f{R) gravity can exhibit "chameleon" behavior: the 
effective mass of the scalar depends on the curvature and, therefore, it can appear as large at 
small scales, making the mediated force short ranged and helping the theory avoid solar system 
constraints; however, the effective mass can become small at large scales, making the force long 
range and allowing the theory to have cosmological effects |741 [TSj [76l I77j . In any case, the 
models have to be very carefully selected in order to be viable and at the same time lead to 
some interesting late time phenomenology (some examples do exist, e.g. |77[ [T8j ) 

The results for Palatini f{R) gravity are much more intriguing. In this case, the known results 
from Brans-Dicke theory cannot be used because the = —3/2 case is exceptional in the sense 
that the scalar is non-dynamical (the standard treatment assumes 7^ —3/2). In [701 171j the 
PPN expansion for ujq = —3/2 Brans-Dicke theory and, consequently. Palatini f{R) gravity 
was rigorously performed. Remarkably, as is pointed out there, and also in [72] without the use 
of the equivalence between theories, whether a model will lead to the correct post-Newtonian 
behavior depends on the density. This had not been encountered before in alternative theories 
of gravity. Just to give a characteristic example, the zero-zero component of the PPN metric is 
given by 

h^S{t,x) = ^-^^ + ^r- + a{n (22) 
r 600 

where Vq and are the background values of the potential and the scalar respectively 
and Mq = (pQ J d^x' p {t, x') /</>. The effective Newton constant G^s and the post-Newtonian 



parameter 7 are defined as 
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7 = 



Mq - My 

Mq + Mv 



(23) 



where My = ^"^0 / d^x' [Vo/cpo - V{(j))/(j)]. The striking term in eq. ([21]) is the last one: ^(r) 
is some algebraic fmiction of the trace of the stress-energy tensor T whose functional form 
depends on the functional form of /. T is in turn algebraically related to the density p and 
so the PPN metric picks up an algebraic dependence on the density! This signals a serious 
shortcoming of Palatini f{R) gravity, non-cumulativity, which will be discussed in the following 
section. Therefore, I refrain from saying more now. 

To summarize, the PPN expansion for f{R) gravity, even though it was a tedious procedure, 
deepened our understanding of post-Newtonian gravity in general and at the same time provided 
important viability constraints. 

2.6. Lesson 6: Non-cumulativity as the root of all evils! 

We just saw that in Palatini f{R) gravity the PPN metric depends algebraically on the density. 
This issue turns out to be related with two other problems: a conflict between Palatini f{R) 
gravity and the standard model of particle physics jTH [80] and the appearance of singularities 
on the surface of spherically symmetric matter configurations [Hll [82l [83] . 

More precisely, it was shown in [79] assuming that matter is a Dirac field and in [80] assuming 
that matter is a scalar field that, when one goes to the local frame where spacetime is flat plus 
second order corrections, non-perturbative corrections and strong couplings appear in the matter 
action even at low energies. Since non-gravitational experiments are held at a local frame, these 
deviations from the standard matter actions contradict the results of such experiments. In 
|81 t I82 t [55] on the other hand, spherically symmetric matter conflguration where considered. It 
was shown that if one assumes a polytropic equation of state and attempts to match any interior 
solution to the unique Schwarzschild-de Sitter exterior, then singularities appear at the surface 
of the matter configuration when the polytropic index F is in the range 3/2 < F < 2. The 
presence of the singularities is essentially insensitive to the choice of the function f{R) |81[ I82j. 
Even though a polytrope is a very idealized form of matter, the importance of this results lies on 
two factors: a) There are physical matter configurations, such as a degenerate non-relativistic 
electron gas or an isentropic monoatomic gas, which can be described accurately by a F = 5/3 
polytrope; such configurations do not have a description in Palatini f{R) gravity, rendering the 
theory incomplete, b) The polytrope is just used as an idealized approximation so that analytic 
results can be easily obtained, and the root of the problem lies deep within the structure of the 
theory. 

Lets us elaborate on this second point. This will also allow us to make the connection with 
the other problems discussed here. As mentioned earlier, the independent connection in Palatini 
f{R) gravity is an auxiliary field, and one can express it algebraically in terms of the metric and 
the matter fields using eqs. ([8]) and ([9]) and eliminate it from the field equations. Indeed, if this 
is done, eqs. ([6]) and ([Tj) reduce to a single equation: 



where G^p is the Einstein tensor of the metric Qfj^p. Note that, since IZ is algebraically related 
to T, f{'TVj and f'(TZ) are essentially function of T. T contains usually first derivatives of the 



G. 




matter fields (or is algebraically related to the density for a perfect fluid). These imply that 
the right hand side (rhs) of eq. ()24|) contains third derivatives of the matter fields (or second 
derivatives of the density). However, this equation is only second order in the metric. Therefore, 
the metric will have some algebraic dependence on the matter fields or the density! Gravity will 
not be cumulative, as in GR or most theories, and the characteristics of the gravitational field 
will depend on the local matter content. 

This explains the behavior of the PPN limit. It also explains the conflict with particle physics 
(the terms on the rhs of eq. (|24p can be considered as coming from strong interactions). An easier 
way to think about this is to recall that in the equivalent Brans-Dicke theory the scalar is non- 
dynamical and algebraically related to the matter but non-minimally coupled to gravity. This 
leads to the strong couplings: when going to the local frame (p cannot be treated perturbatively. 
One has to solve for it in terms of the matter and replace back, introducing the new terms. 
I.e. even though (f) appears not to couple to matter, it enters the matter action though the 
coupling with the metric. 

The singularities found in [8H [82l [83] are a manifestation of the non-cumulativity as well, 
in the sense that, in Palatini f{R) gravity, the metric is bound to inherit any discontinuities 
in the matter and this would lead to a singularity. The simple example of polytropic matter 
configurations is a good demonstration of this, as the range 3/2 < P < 2 for the polytropic index 
is exactly the range for which the density on the surface is finite (and so there is a discontinuity 
since the density is zero in the exterior). However, matter is classically allowed to exhibit 
discontinuities. Remarkably, as it was mentioned, the singularities appear mathematically for 
essentially all functions / (besides GR). There is only one case when one can question their 
physical significance: if they manifest themselves at scales for which the microstructure of matter 
has to be taken into account and the averaging procedure leading to a perfect fluid or a polytropic 
description is questionable. In this sense requiring the absence of such singularities provides a 
very strong constraint, ruling out all models with late time cosmological phenomenology, but 
even models with ultraviolet corrections [81^ I82j . Only models with Planck suppressed correction 
seems to survive |84j . 

3. Conclusions 

As is probably obvious by now, proposing and studying gravity theories other than GR is not an 
easy task. One can easily end up causing more problems that those the proposed modification 
of gravity is attempting to solve. Many of the shortcomings of f{R) gravity become much more 
obvious after field redefinition and representation changes. The fact that it took some time to 
discover them seems to indicate that habit and affinity to specific formalisms and representation 
of a theory seem to have played a crucial role in the misunderstandings and the difficulties of 
f{R) gravity. In any case, f{R) gravity has been a very successful toy theory. It seems that 
we may have to go beyond simple f{R) models if we want to have a chance of explaining some 
of the current riddles of physics with a deviation from GR. However, our experience with f{R) 
gravity indicates that pursuing such an idea might be fruitful: possibly in providing a solution 
to some specific open problem, but most certainly in providing insight into the most puzzling 
interaction, gravity. And this is an important lesson on it own. 
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